Comparison of Propagator Decomposition in Seismic Imaging by Wavelets, Wavelet-packets, and Local Harmonics
نویسندگان
چکیده
Kirchhoo migration operator is a highly oscillatory integral operator. In our primary work 1] (Wu and Yang, 1997), it has been shown that the matrix representation of Kirchhoo migration operator for homogeneous background in space-frequency domain is a dense matrix, while the compressed operator in beamlet-frequency domain, which is the wavelet decomposition of the Kirchhoo migration operator, is a highly sparse matrix. Using the compressed matrix for imaging (beamlet migration), we can retain the wide eeective aperture of a full-aperture operator, and hence achieve higher resolution and image quality with reduced computational cost. However, as well known, wavelets work best for zero-frequency stationary signals. But, for the Kirchhoo migration operator, it is both space-varying in the near-eld region and high-frequency stationary in the far-eld zone. Therefore, wavelets are not very eecient for this kind of operator. In this research, we rst summarize the results of maximum sparsity adapted wavelet-packet transform (MSAWPT) for the decomposition and compression of Kirchhoo migrator 2] (Wang and Wu, 1998), and then further study the decomposition and compression of Kirchhoo migration operator by local harmonics (i.e., local cosines/sines). It was found in 2] (Wang and Wu, 1998) that the MSAWPT can generate a more eecient representation for the imaging operator than the standard discrete wavelet transform (DWT) and the compression capability of MSAWPT is much greater than that of DWT. In this paper, we also observed that for low frequency operator, the compression capability of uniform local cosine bases is equivalent to that of standard wavelets and is weaker than that of adapted wavelet-packets; while, for high frequency operator, uniform local cosine bases are more powerful than both the standard wavelets and adapted wavelet-packets. Furthermore, for local cosine transform, a good parameter setting (i.e., type and smoothness of the bell function, edge extension, overlapping radius, folding style, and window size) can generate a higher compression ratio.
منابع مشابه
A021 Compression Denoising: Using Seismic Compression for Uncoherent Noise Removal
Several authors, e. g. Vermeer et al. [1996], have recently recognized that wavelets may not be the best fit for seismic data since they present large-scale oscillations. They have therefore investigated for instance wavelet packets, local cosine or Gabor transforms. We focus here on filter banks (FB), which can be regarded as a multichannel generalization of wavelets [Duval and Røsten, 2000]. ...
متن کاملDiffusion wavelet packets
Diffusion wavelets can be constructed on manifolds, graphs and allow an efficient multiscale representation of powers of the diffusion operator that generates them. In many applications it is necessary to have time–frequency bases that are more versatile than wavelets, for example for the analysis, denoising and compression of a signal. In the Euclidean setting, wavelet packets have been very s...
متن کاملA combined Wavelet- Artificial Neural Network model and its application to the prediction of groundwater level fluctuations
Accurate groundwater level modeling and forecasting contribute to civil projects, land use, citys planning and water resources management. Combined Wavelet-Artificial Neural Network (WANN) model has been widely used in recent years to forecast hydrological and hydrogeological phenomena. This study investigates the sensitivity of the pre-processing to the wavelet type and decomposition level in ...
متن کاملA Study on B-Spline Wavelets and Wavelet Packets
In this paper, we discuss the B-spline wavelets introduced by Chui and Wang in [1]. The definition for B-spline wavelet packets is proposed along with the corresponding dual wavelet packets. The properties of B-spline wavelet packets are also investigated.
متن کاملFault Strike Detection Using Satellite Gravity Data Decomposition by Discrete Wavelets: A Case Study from Iran
Estimating the gravity anomaly causative bodies boundary can facilitate the gravity field interpretation. In this paper, 2D discrete wavelet transform (DWT) is employed as a method to delineate the boundary of the gravity anomaly sources. Hence, the GRACE’ satellite gravity data is decomposed using DWT. DWT decomposites a single approximation coefficients into four distinct components: the appr...
متن کامل